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ABSTRACT 

G73.9+0.9 has been elassified as a probable shell-type supernova remnant (SNR), 
although it has also been suggested that this objeet eould be a pulsar wind nebula 
(PWN). Here, a broadband model of the non-thermal emission of G73.9-I-0.9 from 
radio to gamma-rays is presented. The model ineludes a new gamma-ray observation 
obtained with analysis of 7 years of data from the Eermi LAT teleseope. Above 200 
MeV the source is detected with a significance of I3a and the spectrum of the radiation 
is best described by a power law with an index of ~2.5. The leptonic mechanisms are 
hard to reconcile with the measured radio and gamma-ray SED. A PWN origin for the 
high-energy emission is also not very likely, due to the lack of detection of pulsars and 
of X-ray emission in the region, as well as from the shape of the gamma-ray spectrum. 

Given the possibility that the object is interacting with molecular clouds, a hadronic 
origin of the high-energy emission is more likely, and the spectral properties of the 
cosmic rays responsible for this radiation are derived. 

Subject headings: acceleration of particles-radiation mechanisms: non-thermal-ISM: 
individual objects: G73.9-I-0.9-ISM: supernova remnants. 


1. Introduction 

In the last decade, an increasing number of supernova remnants (SNRs) have been detected 
at gamma-ray energies. Establishing this type of objects as high-energy emitters is important as 
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photon signatures could help solve the long-standing problem on the origin of Galactic cosmic 
rays, charged particles with energies up to 10^^ eV. These particles are thought to experience ac¬ 
celeration in the shocks of SNRs which could transfer up to 10% of the explosion kinetic energy 


Bell 1978: Blandford & Eichler 1987) 

. Observations of synchrotron radiation from SNRs, from 

radio to X-rays ( 

Gotthelf et al. 2001: 

Berezhko et al. 2002: 

Hwang et al. 2002: 

Rho et al. 2002: 

Lone et al. 2003: 

Vink & Lamina 20C 

)3: 

Berezhko & Volk 2004). reveal the existence of high- 


energy electrons in these sites. 

Gamma-rays may be produced by high-energy particles mainly through the mechanisms of 
inverse Compton (IC) up-scattering of low energy ambient photons, non-thermal bremsstrahlung 
radiation and the decay of neutral vr—mesons produced in inelastic collisions between accelerated 
ions and ambient nuclei. The hadronic mechanism may produce a gamma-ray spectrum which 
can be hard to distinguish from leptonic emission. Some SNRs are bright gamma-ray sources and 


act with these clouds ( 

Abdo et al. 2010a: 

Abdo et al. 2010b: 

Abdo et al. 2010e: 

Abdo et al. 20101: 

Abdo et al. 2009: 

Ackermann et al. 2013 

: Aiello et al. 2012 

: Ixing et al. 2014 

Xing et al. 2015: 


Liu et al. 


m. These observations support the idea that cosmic rays are produced in SNRs. 
The high-density material in MCs could provide abundant targets for gamma-ray production in 
a hadronic scenario, making SNR-MC interactions ideal targets to constrain the spectrum and 
energetics of cosmic rays. On the other hand, the origin of gamma-rays in younger remnants re- 
mains controversial. While the emission from Tycho’s SNR or Gas A might be of hadronic origin 
( Abdo et al. 2010d ^ Arava 201^: Araya 2011), others show gamma-ray spectra that is charact eris- 
tic of IC emission ( Ellison et al. 2010 : Yuan et al. 2011 : iTanaka et al. 2011 : Yuan et al. 2014 ). 


At high energies, SNRs have been observed by the Fermi LAT, a pair conversion telescope 
with a field of vie\y of 2.4 sr sensitive to gamma-rays from ~ 20 MeV to more than 300 GeV 
( Atwood et al. 2009 ). Here, a detailed analysis of LAT data in the region of the SNR G73.9-I-0.9 
and broad-band modeling are presented. This object is located near the complex Cygnus region 
and it is seen in radio images as a compact nebulosi ty with a size of ^ 27'. It is classified as an SNR 
based on its polarization and nonthermal spectrum (iReich et al. 1986). Its measured flux density is 
9 Jy at 1 GHz and the photon spectral index is (F^, oc u~°‘) a = 0.23 ( Kothes et al. 2006 1. The type 
of remnant is uncertain, the possibilities include a shell type remnant or a pulsar wind nebula. No 
pulsars were found in a search by Gorham et al. (119961) wit h a sen sitivity of 0.3 mJy at 430 MHz 
and 1400 MHz, and another pulsar search by Lorimer et al. (Il998h also resulted in no detection. 


The distance to the source is relatively unknown. Lozinskaya et al. (119931) derive a distance 
of 1 — 2 kpc based on optica l observations. B ased o n the S — 79 relation, a distance of ~ 6 kpc was 
determined by Leahy (119891) . Mavromatakis (120031) conclude that a distance of 2 kpc is compatible 
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with t heir op tical observations of the source. Lorimer et al. ( 1998 ) cite the study by Lozinskaya 
et al. 1 1993 b and use a distance of 1.3 kpc. Pineault et al. (Il996h find an interesting HI region 
that shows a correlation with the remnant’s morphology at a kinematic distance of 10 kpc, but this 
distance is probably less likely for G73.9+0.9. A distance closer to more recent estimates of 1.8 
kpc is adopted here. 


From their observations, Lozinskaya et al. (1 1993b also estimate a shock expansion velocity of 
200 — 300 km/s and age of 11 — 12 kyr for G73.9+0.9, for an ambient density of 10 cm“^. They 
compare these properties to those of the well-known SNR IC 443 and find similarities. 


A number of compact radio sources are seen near the remnant, most of w hich have spectra tha t 
are compatible with HII regions or dust pressumably heated by the shock (iPineault et al. 1990b . 

J = 1 — 0 observations show molecular clouds with spatial correlat ions w ith the radio 
emission from the SNR at small CO velocities (+1 to +3 km/s, Jeong et al. (l2012b l to the north 
of the radio peak, although no direct evidence for an interaction has been reported. Other studies 
of the ambient near the SNR have concluded that there is an HI shell around it that are likely 
associated w ith it, and have also detected molecular cloud emission to the northwest of G73.9+0.9 
( Sitnik 2010 ). 


Mavromatakis (I2003h derive column densities that imply ISM densities of 0.5 cm“^ at the 
distance of 2 kpc and conclude, based on the properties of the optical emission, that an extended 
radio blob seen to the north of the SNR is probably unrelated to it, but pe rhaps associated to a 
bipolar HII region. This had also been concluded by Pineault et al. (Il996b after comparing the 
infrared properties of the SNR and the northern component, which ar e very different, as is the 
morphology of the HI gas in both regions. According to Mavromatakis (l2003b . the low ionization 
optical images further indicate an inhomogenous and patchy intestellar medium. 


X-ray observations fro m the ROSAT All-Sky survey of the region show no especially promi¬ 
nent soft emission. Leahy (Il989b analysed Einstein observations and derived a 2a upper limit on 
the 0.2 — 4 keV flux of 3.0 x 10“^^ erg cm“^ s“^. There is an archival XMM-Newton observation 
of the remnant that was proposed to search for a pulsar wind nebula but no obvious emission was 
found (private communication). At very-high e nergies (VHE), VERITAS set a 99% upper limit 
for the source above 300 GeV ( Theiling 2009 ) of 2.33 x 10“^^ photons cm“^ s“^ with ~ 600 
minutes of data (assuming a power-law spectrum, dN/d E oc E~^, with index T = 2.5). The 
radio fluxes use d in this work are taken from the literature dReich et al. 1986 : Pineault et al. 1990 : 


Sun et al. 2011 ) 


In this work, the gamma-ray observation is presented in Section and its properties explored. 
The spectral energy distribution (SED) is modelled in Section [3l Different possibilities for the 
origin of the emission are studied. Finally, the conclusions are given in Section HI 
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2. LAT Data 


LAT data gathered from the beginning of the mission, 2008 August, to 2015 July were 
analysed with the most reeent software, SCIENCETOOLS version vlOrOpil released 2015, June 
24, with the latest reproeessed “Pass 8” photon and spaeeeraft data and the instrument response 
functions P8R2_SOURCE_V6P. The Pass 8 is an effort to improve the reconstruction algorithms 
applied to the data and includes an overall improvement in the understanding of the detector 
( Atwood et al. 2013 ). 

Photons with energies from 200 MeV to 100 GeV are selected within a 14° x 14° square region 
of interest (ROI) centered at the position of G73.9+0.9, (J2000) a = 303.54, 6 = 36.20. Below 200 
MeV the steeply falling effective area and broadening of the energy dispersion and PSF present 
significant challenges in the analysis, and above 100 GeV the drop in the number of photons 
becomes an issue. 


The data selection is limited to zenith angles less than 90° to remove photons from the Earth 
limb and the event class is cut to keep only the SOURCE class events as recommended. Time 
intervals are selected when the EAT was in science operations mode and the data quality was 
good. The data is binned in 30 logarithmically spaced bins in energy and a spatial binning of 0.°1 
per pixel is used. 


The anal ysis i ncludes all sources listed in the latest Fermi EAT catalog (hereafter 3FGE, 
Acero et al. (l2015h f. as well as the standard diffuse Galactic emission model and the isotropic 
model that accounts for the extragalactic gamma-ray background and misclassified cosmic rays, 
both of which are provided with the analysis tools. The normalizations of these spatial templates 
are left free during the likelihood fit, as well as the normalizations of the cataloged sources located 
less than 6° from the center of the ROI, to account for the broad PSF at low energies. The spectral 
parameters of all the sources are kept fixed during the fit. 


The properties of EAT sources are obtained by a maximum likelihood analysis q which esti¬ 
mates the probability of obtaining the data given a source model and fits the source parameters to 
maximize this probability. The unidentified source candidate 3FGE J2014.4+3606, located about 
0°.l from the center of the ROI, is reported in the 3FGE catalog with a significance of 4.5a. This 
source candidate is not included in the model in order to accurately study the morphology of the 
residuals and the significances of new sources seen with the larger data set used here. 


The bright EAT pulsar 3FGE J2021.1+3651 is found l.°6 from the center of the ROI. This 


' See http;//fermi.gsfc.nasa.gov/ssc 

^See http;//fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/binnedJikelihood_tutorial.html 
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source shows a spectrum with a cut off at 3 GeV. Other nearby sources include the FSRQ 3FGL 
J2015.6+3709 (l.°l from the ROI center, also known as VER J2016+371) and 3FGL J2017.9+3627 
(located ~ 1° from the SNR), believed to be associated to the VHE source MGRO J2019+37. 
These sources have EAT spectra that are described by a logParabola. In order to investigate the 
effect on the residual emission near these sources, fits are done in different energy intervals keeping 
these source parameters free, and the resulting residuals maps obtained after subtracting the best- 
fit model to the data are compared to those obtained when their spectral parameters are fixed to 
the values reported in 3EGE. The results are very similar in both cases. The spectral analysis in 
Section llTl was also repeated keeping the spectral shape of nearby sources free without significant 
differences in the results. The results of the analysis reported here was then obtained with spectral 
parameters of all sources fixed to the values found in the EAT catalog. 


In order to remove the emission of nearby gamma-ray sources as well as the diffuse back¬ 
grounds, the tool gttsmap is used above 3 GeV in a 1.2° x 1.2° region around the center of the ROI. 
Eor each pixel of the map the test statistic (TS) is evaluated. The TS is defined as —2 log(Lo/T), 
where Lq and L are the maximum likelihood values for the null hypothes is and for a model in clud- 
ing an additional point source at the location of that pixel, respec tively (IMattox et al. 19961) . The 
detection significance of a source is given by \/TS (however see IProtassov et al. 20021 for some 
caveats) and thus a TS value of 25 is considered the threshold for source detection. 


Pig. [Ushows the resulting significance map around G73.9-I-0.9 as well as the locations of other 
objects in the region. The PSP of the EAT suggests that the excess found here is consistent with a 
point source, as it would be expected for an object with the extension of SNR G73.9-I-0.9, or any 
structure within it. In order to find the position of the new source the tool gtfindsrc is applied to the 
same data set above 3 GeV, the resulting best-fit coordinates are (J2000) a = 303.337, 6 = 36.226 
with a 68% confidence level positional error radius of 0°.05. As an additional test, residuals maps 
are also obtained at lower energies and the excess observed seems consistent with a point source 
at this position, indicating that there are no further sources in the region having a softer spectrum 
that would not have shown up in the TS map above 3 GeV. 


The map shows a separation of ~ 0°.25 between the new best-fit position and the position of 
the source candidate reported in the 3PGE catalog, obtained with a smaller data set but also with 
different analysis software versions and diffuse background models. The significance of the new 
source reported here is also different (~ 6a above 3 GeV, but much higher for the whole energy 
range, see below). A satisfactory substraction of the emission is also seen in the residuals map 
once this source is added to the model. 
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2.1. Spectrum of the source 


Once the new souree position and morphology are known, its speetrum is studied above 200 
MeV with the likelihood teehnique. Different speetral shapes are probed in order to eompare the 
resulting TS values: simple power-law, broken power-law, logParabola, and power-laws with ex¬ 
ponential and ‘super’ exponential eut offs. It is found that models that require more free parameters 
with respeet to the simple power-law result in similar TS values, whieh are (in the same order): 
177, 184, 180, 182 and 181. Thus, the fits with all these funetions eorrespond to an improvement 
at the ~ 2cr level with respeet to the simple power-law, whieh is not very signifieant. The best-fit 
speetral shape adopted in the rest of the analysis is a simple power-law. The values of the integrated 
souree flux for these different speetral shape fits are very similar. 


Above 200 MeV the speetral index of the souree is (2.51 ± O.OTstat) with an integrated flux 
above this energy of (1.87 ± 0.21stat) x 10“® em“^s“^, and the overall signifieanee of the souree is 
quite high, IScr. The eorresponding luminosity above 1 GeV is (2.7 ± O.Sstai) x 10^^ g erg s“^, 
where di g is the distanee in units of 1.8 kpe. This value is very similar to th e luminosity of other 


dim gamma - ray remnant s sueh as the Cygnus Loop, G296.5-1-10.0 and HB9 (iKatagiri et al. 2011 


Araya 20131 : lAraya 2014h . The souree TS drops to 3 for events above ~ 6 GeV. 


In order to eonstruet a gamma-ray SED the data is binned in 9 logarithmieally-spaeed intervals 
from 0.2 MeV to 100 GeV, and a likelihood fit is applied within eaeh keeping the speetral index of 
the souree fixed to the global value. Whenever the TS of the souree falls below 9 a 95% eonfidenee 
level upper limit is derived in the eorresponding interval. Several systematie uneertainties are 
taken into aeeount in the SED fluxes for eaeh bin. Aeeording to the analysis of Pass 8 EAT data, 
the uneertainty on the effeetive area is about 5% and the systematie uneertainty on the PSE is 
always less than 5% in the energy range used here. These uneertainties are added in quadrature to 
the flux statistieal errors in the SED. 


2.2. Variability analysis 


Extragalaetie gamma-ray sourees are known to be highly variable. In order to seareh for 
variability, the data is binned in ten ~ 8-month bins and the likelihood analysis is applied to 
estimate the flux of the souree seen in the direetion of G73.9-I-0.9. The normalizations of the other 
sourees in the model are kept fixed to the values found for the fit to the 7-yr data set, as well as 
the speetral index of G73.9-I-0.9. Eig. [2] shows the resulting fluxes for eaeh bin. A fit to a eonstant 
value gives an average flux of (1.64 ± 0.14) x 10“® om“^s“^ (la error), and a redueed of 1-9. 


In ord er to quantify the probability of ehanges in the flux with time, a variability index is 
ealeulated (lAbdo et al. 201 Pel) and the result is 17.2 eonsidering statistieal errors on the data only. 
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Varibility is considered probable when this index exeeeds the threshold of 21.7, eorresponding to 
99% eonfidenee in a distribution with 9 degrees of freedom. Adding systematie errors would 
result in a lower index. It ean be said that the souree shows no signifieant long term variability. 


3. Discussion 

Fermi LAT data eonfirm the existenee of signifieant gamma-ray emission from an unresolved 
souree at the loeation of G73.9-I-0.9, an objeet whose nature is still debated. The speetrum and 
properties of the souree are now studied in the eontext of multi-wavlength obervations eonsidering 
several seenarios. 


3.1. A pulsar or pulsar wind nebula 


Pulsars are an important elass of gamma-ray sources. No radio pulsars have been found in 
the region of the SNR. Most pulsars have gamma-ray s peetra that ean be fit with a flat power-law 
that euts off ex ponentially at energies of a few GeV (lAbdo et al. 20131: lAekermann et al. 20121: 
Lee et al. 2012h . while the speetrum of the gamma-ray souree deteeted here has been shown to be 
deseribed by a simple power-law. Some pulsars reeeive a “kiek” during the supernova explosion. 
The angular separation from the position of the gamma-ray souree and the peak of the radio shell 
is ~ 0.°14. Even if the remnant is assumed to be relatively young (age ~ 1 kyr), which is unlikely 
given the laek of detection of X-ray emission, the corresponding pulsar transverse veloeity is 4500 
km/s for a distance of 1.8 kpe, whieh is high but still possible; however, no pulsars have been found 
in the region. 

Pulsars transfer part of their rotational energy to a relativistie magnetized wind (mostly eleetron 
positron pairs). When this wind slows down abruptly it ereates a termination shoek where par- 
tiele aeeeleration oeeurs, resulting in a pulsar wind nebula (PWN). PWNe are seen from radio 
to gamma-rays. Their non-thermal emission is dominated by synehrotron an d inverse Compton 
from high-energy leptons, although hadronic processes may also contribute dCheng et al. 19901: 
Amato et al. 2003 ). 


G73.9-I-0.9 has several radio properties that are similar to those of PWNe such as a low spectral 
index. However, PWNe deteeted by the LAT typieally have young, bright pulsars, and are also seen 
at TeV energies (lAbdo et al. 201 In . Hard X-ray and gamma-ray observations of PWNe have the 
potential to eonstrain, among other things, the speetral break and maximum energy of synehrotron- 
emitting partieles. The radio speetrum of a PWN is hard {Fy oc z/“", a = 0 — 0.3) and softens 
at higher energies, and in the X-ray band a > 1. The synehrotron speetrum of the Crab Nebula, 
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for example, may be described as power-laws with breaks around 10^^ Hz, 10^® Hz and 100 keV 
and only the peak at VHE is attributed to IC emission. PWNe usually have relatively high radio 
polarization fractions (from 30% to 50% at 1 GHz), while for G 73.9-r0.9 it is closer to the value 
found for SNRs, around 1% to 5% at ~ 5 GHz ( Sun et al. 2011 ). 


PWNe can also have a wide range of nebular magnetic fields (~ 5 /iG to > 1 mG, e.g., 
Reynolds et ah 20121) and show a variety of break frequencies in their radio spectra. Fig. [3] shows 
the SED from radio to gamma-rays with two possible models that include synchrotron emission in a 
magnetic field of 4 /iG and IC scattering of CMB photons (this low field value is only a lower limit, 
as the IC emission from electrons up-scattering other ambi ent photons is not considered). The ra¬ 
dio fluxes corre spond to integrated emission of G73.9-1-0.9 (IReich et al. 19861 : iPineault et al. 1990l: 
Sun et al. 20111) . The Einstein UE at X-ray energies implies that the GeV emission cannot be pro¬ 
duced by synchrotron emission from the same radio-emitting particles. The models are shown for 
broken power-law particle distributions with different indices and break and maximum electron 
energies: 1.46, 4.0, 204 GeV and 1 TeV, and 1.46, 2.2, 50 GeV and 600 GeV, respectively. A 
lepton energy density of ~ 10“^° erg cm“^ is used in both cases, which is much higher than the 
resulting magnetic field energy density (6.4 x 10“^^ erg cm“^). The total particle energy content is 
~ 4 X 10"^® erg and the models assume a source distance of 1.8 kpc. Contributions from synchrotron 
self-Compton emission are negligible. 

No known EAT PWNe show a peak of the IC emission in the EAT energy range. In a search 
for EAT PWNe, the confirmed PWNe and the candidates were always ass ociated with young and 
powerful pulsars with spin-down power between 10^® and 10^® erg s“^ (lAcero et al. 20131) . al¬ 
though the search was biased as the objects studied were previously selected for having ass ociated 
VHE emission. In a systematic exploration of models of Crab-like PWNe, Torres et al. (120131) 
studied a variety of parameters affecting the SED of these objects, including the age. From their 
work, no case is seen where a source would show a very low X-ray flux (< 10“^^ erg cm“^ s“^) 
and a negative SED slope at GeV energies as found here for G73.9-I-0.9, at least up to a source age 
of 9 kyr, similar to the possible age of this object. 


3.2. An AGN? 


Active Galactive Nuclei (AGN), and blazars in particular, are an important class of gamma-ray 
sources. 3FGE includes 1745 sources associated with AGN (58% of all 3FGE sources), and 98% 
of them are blazars or candidate blazars, the rest being radio galaxies, Seyfert galaxies and other 
types of AGN. Most of the gamma-ray emittin g AGN are bright so urces of radio emission. No 
blazars are found in the region near G 73.9-r0.9 (IMassaro et al. 20091) and no other AGN types are 
found in the Million Quasars Catalog (iFlesch 2015h (version 4.5, May 2015) near the gamma-ray 
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source. The high-energy emission is also found to be consistent with steady emission in ~yearly 
time scales (Fig. |^. 


3.3. SNR-MC interaction 


There is no direct evidence of an interaction of G73.9-I-0.9 with a MC and more studies are 
necessary to confirm or reject this sc enario, but the morphological properties of the ambient clouds, 
as observed by Jeong et al. (l2012h . point to this possibility. In particular, there is an interesting 
compact MC within the best-fit position of the gamma-ray emission, as can be seen in Fig. [TJ 


Fig. m shows an SED model for which the dominant contribution to the gamma-ray emis sion 
is hadronic, calculated with a recent parametrization of the cross section (iKafexhiu et al. 20141) for 
these interactions. The cosmic ray distribution used to calculate the emission is a simple power-law 
in momentum of the form The use of this proton distribution is motivated by the shape of 

the gamma-ray spectrum. Hadronic gamma-ray emission is expected to increase rapidly around 

photon energies of 100 MeV and then follow the particle distribution above ~ 1 GeV. The model is 

/ \ -1 


calculated for a total proton energy of 10^° V icm~^ j target proton density. 

Since MCs can have a wide range of densities (including ?7,p >> 1 cm“^), the required total particle 
energy is realistic. The proton power-law spectrum in the model cannot extend beyond a particle 
energy of 15 — 20 GeV, based on the EAT upper limits. 


Even though the SED obtained with this model reproduces the observed SED, the choice of 
a broken power-law for the particle distribution can also succesfully explain the data. Gamma- 
ray spectra observed in systems with SNR-MC interaction s are sometimes nrqdeled with a bro¬ 
ken power-law ( e.g., Abdo et al. 2009 : Abdo et al. 2010e : Abdoetal. 20101 : Ohiraet al. 2011 : 


Tang et al. 201 ih . In this scenario, a proton distribution with a break momentum of ~3 GeV c 


-1 


and indices of 1.5 below (as found for the leptonic component) and 3 above this break, with a total 

cosmic ray energy of 0.8 x 10®° ( icm~° ) reproduce the data as well. These parameters, 

however, are not very well constrained, as also lower break momenta and lower indices above the 
break (i.e., harder proton distributions) are consistent with the data. According to previous work, 
the break i nomentum depends on th e magnetic field, ion density and the frequency of ion-neutral 
collisions (IMalkov & Sagdeev 201 ih . but these parameters are unknown. 


Eor the model shown in Eig. HI the leptonic component is obtained with a simple power-law 
for the electron distribution (index 1.46 and maximum electron energy of 50 GeV) with a total 
energy in leptons of 1.5 x 10^^ erg and a magnetic field of 30 /iG. If the field is closer to the 
average Galactic value, ~3 /iG, the corresponding total lepton energy would be 2.8 x lO"^® erg. 
Such low fields would still be sufficiently capable of confining protons with energies of tens of 
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GeV producing the gamma-rays. 


3.4. Other leptonic mechanisms: hremsstrahlung and IC emission 

It is difficult to explain the LAT fluxes below 400 MeV with IC-CMB emission from electrons 
following a power-law distribution with index 1.46. The magnetic field used is 12 fxG and the 
maximum electron energy is 500 GeV. The synchrotron SED emission would peak around infrared 
frequencies (~ 3 x 10^^ Hz) at the level of ~ 9 x 10“^^ erg cm“^ s“^. This scenario can be seen 
in Fig. [51 

The relatively hard electron spectrum in this source makes a bremsstrahlung-dominated sce¬ 
nario for the gamma-rays very unlikely. Fig. [5] also shows this case for an ambient target density 
Up = 1 cm“^. The predicted gamma-ray flux could become somewhat similar to the level mea¬ 
sured only when the maximum electron energy is well below ~ 30 GeV, a value for which the 
resulting radio spectrum starts to deviate from the data. The hremsstrahlung scenario can be safely 
discarded for the gamma-rays in G73.9-I-0.9, as it is unable to explain the gamma-ray data. 


4. Conclusions 

Gamma-ray emission has been found in the region of G73.9-I-0.9 with cumulative observations 
by the FAT onboard Fermi. The spectrum of the emission can be best described with a power-law 
with index 2.5 in the interval where the signal is significant, between 200 MeV and ~6 GeV. No 
other known possible gamma-ray emitter is seen at or near the location of the point-like source. 
Being in the region where there is a possible SNR, it is natural to explain the origin of the gamma- 
rays as due to interactions of high-energy particles with the ambient gas or radiation fields, but 
there are other possibilities. 

Several scenarios were considered to explain the origin of the emission. There are no known 
pulsars in the region. The gamma-ray spectrum of pulsars usually shows a cutoff or curvature 
around a few GeV, which is not seen for this source. Also, there are no AGNs that are known in 
the direction of G73.9-I-0.9, and the flux from the source was shown to be consistent with steady 
emission in a time scale of years. The lack of detection of X-ray emission and the gamma-ray 
spectral shape are in conflict with the PWN scenario for the broadband SED, according to models 
for PWN with ages up to 9 kyr. 

Due to all these constraints, it is more likely that the gamma-rays originate from interactions of 
high-energy particles accelerated in the SNR. The shape of the SED implies that the IC scenario is 
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not very probable. On the other hand, the bremsstrahlung seenario is ineonsistent with observations 
and ean be disearded. 


From the spatial eorrelation between the SNR and CO emission from a moleeular eloud, 
whose presenee would be ideal to enhance the gamma-ray luminosity from hadronic interactions, 
this possibility seems more favorable. CO emission is indeed seen within the best-fit position 
of the gamma-ray source. The GeV SED of this source is also consistent with a hadronic sce¬ 
nario and, within the uncertainties, the particle population responsible for the emission can be 
described by either a simple power-law (index of 2.5) or a broken power-law in momentum. Bro¬ 
ken power-law distribut ions have been observed for the particles when the sho ck of a SNR in¬ 
teracts with a MC (e.g., Abdo et al. 2009 : Abdo et al. 2010e : Abdo et al. 20101 ). The steepening 
in the distribution might be due to poor particle confine ment in regions with low ionization (e.g.. 


Ptuskin & Zirakashvili 20031: IMalkov & Sagdeev 201 ih . If, as is the case for the leptonic popu¬ 
lation, a hadronic index of ~1.5 is used for the first component of the power-laws, the required 
index is 3 above a particle momentum break of 3 GeV/c, however, these parameters are not well 
constrained by the data. 
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Fig. 1.— Significance map of the gamma-ray emission around G73.9-I-0.9 output from gttsmap 
above 3 GeV, in units of test statistie (TS, see text). The map has a dimension of 1.2°x 1.2° and 
a resolution of 0°.06 per pixel. The magenta eontours mark the non-thermal radio emission of 
the SNR for five eq ually-spaced levels between 0 and 0.3 Jy/beam, taken from the 4850 MHz 
Green Bank survey (ICondon et al. 19941) . The center of the cyan circle indieates the best-fit po¬ 
sition (and its radius the positional uneertainty) of the gamma-ray souree, the white eross indi¬ 
eates the position of the previous 3FGL souree eandidate and black diamonds denote the positions 
of radio knots whieh are likely HII regions. The blue eircle marks the loeation of the peak of 
emissi on from a molecular cloud, obtained from J = 1 — 2 observations by Jeong et al. 
(1201 211 ■ No other 3FGL souree is found within this map and no known blazars are seen in it ei¬ 
ther (IMassaro et al. 200911 . The yellow circle shows the approximate size of the LAT PSF for 68% 
containment of FRONT-i-BACK events. 
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Fig. 2.— Lightcurve of the source seen at the location of G73.9+0.9. The fluxes are calculated for 
10 time bins in the energy range 0.2 — 100 GeV and show only statistical errors. The solid line 
corres lue. 




Fig. 3.— SED of G73.9+0.9 and two simple one-zone leptonic models in the PWN scenario 
(synchrotron plus IC-CMB). The gamma-ray points and upper limits are from the LAT observation 
in this work. The particle distribution in both models is a broken power-law with indices and break 
particle energies of: 1.46, 4.0, 204 GeV (solid lines); and 1.46, 2.2, 50 GeV (dashed lines). 
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Fig. 4.— Hadronic model for the SED of G73.9+0.9. The solid line eorresponds to synchrotron 
emission (magnetic field of 30 /iG, maximum electron energy 50 GeV), the dashed line is the 
gamma-ray flux from eosmie-ray interactions (ambient density Up = 1 cm“^, total eosmic ray 
energy 10^° erg) and the dotted line IC scattering of CMB photons by the radio-emitting electrons. 
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Fig. 5.— Leptonic models for the SED of G73.9+0.9. The solid lines eorrespond to synehrotron 
and IC-CMB in an IC-dominated seenario (a magnetie field of 12 fiG and maximum eleetron 
energy of 500 GeV are used). The dashed lines are synehrotron emission (maximum eleetron 
energy 30 GeV, magnetie field 3 fxG) and non-thermal bremsstrahlung (target ambient density 
Up = 1 cm“^) in the seenario where this latter emission dominates the gamma-ray flux (and for 
whieh the dotted line shows the corresponding IC-CMB flux level). 




